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Abstract 



An optimal control problem associated with the dynamics of the orientation of a bipolar 
molecule in the plane can be understood by means of tools in differential geometry. For first time 
in the literature fc-symplectic formalism is used to provide the optimal control problems associated 
to some families of partial differential equations with a geometric formulation. A parallel between 
the classic formalism of optimal control theory with ordinary differential equations and the one 
with particular families of partial differential equations is established. This description allows us 
to state and prove Pontryagin's Maximum Principle on /c-symplectic formalism. We also consider 
the unified Skinner-Rusk formalism for optimal control problems governed by an implicit partial 
differential equation. 

1 Introduction 

Boscain et al. in [6] study the controllability of the equation 



which models the rotation motion of a bipolar rigid molecule confined to a plane with two control 
electric fields. See references in the paper for more details about the problem origin and its interests. 

The study of controllability in does not approach the problem of existence and construction of 
suitable controls for governing the position of the molecule. The controls are obtained depending on 
the purpose to be achieved. For instance, if the controls are related with the energy needed to take 
the molecule to a particular position or to track a molecule path in the configuration space, then we 
might be interested in minimizing the energy consumption. In other words we have associated an 
optimal control problem to the above control equation. Note that control equation is a particular 
kind of second-order and control-linear partial differential equation. 

In this paper, a geometric approach is considered to deal with optimal control problems for 
particular families of control partial differential equations, similar to the above mentioned example. 

Following the same lines as in our previous paper on Pontryagin's Maximum Principle [3], PMP, 
we extend the geometric method to optimal control problems with some control partial differential 
equations. As in the classical PMP, to succeed in this extension is necessary to extend the control 
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system in a proper way and to optimize a suitable function, the Pontryagin Hamiltonian, instead of 
the functional. 

Before entering into the details, let us provide some historical background on optimal control 
theory. L. S. Pontryagin talked publicly for the first time about the Maximum Principle in 1958, in 
the International Congress of Mathematicians that was held in Edinburgh, Scotland. This Principle 
was developed by a research group on automatic control created by Pontryagin in the fifties. He 
was engaged in applied mathematics by his friend A. Andronov and because scientists in the Steklov 
Mathematical Institute were asked to carry out applied research, especially in the field of aircraft 
dynamics. 

At the same time, in the regular seminars on automatic control in the Institute of Automatics 
and Telemechanics, A. Feldbaum introduced Pontryagin and his colleagues to the time-optimization 
problem. This allowed them to study how to find the best way of piloting an aircraft in order to 
defeat a zenith fire point in the shortest time as a time-optimization problem. 

Since the equations for modelling the aircraft's problem are nonlinear and the control of the 
rear end of the aircraft runs over a bounded subset, it was necessary to reformulate the calculus 
of variations known at that time. Taking into account ideas suggested by E. J. McShane in [14] . 
Pontryagin and his collaborators managed to state and prove the Maximum Principle, which was 
published in Russian in 1961 and translated into English [16] the following year. See [5] for more 
historical remarks. 

Although geometric control theory has been studied since the sixties, geometric optimal control 
theory started to be developed in the nineties jlH I19j . However, no geometric description has 
been made for optimal control problems governed by partial differential equations [10\ \T8\ [20] . In 
this paper we extend for first time the geometric description of optimal control problems to those 
governed by some classes of partial differential equations in order to solve optimal control problems 
that can be found in the physical world, as mentioned above. 

With this purpose in mind, the natural geometric background to use is /c-symplectic formalism 
(Giinther standard polysymplectic), which is a generalization of the symplectic formalism in classical 
mechanics. The fc-symplectic formalism makes possible to geometrically interpret some problems 
such as the vibrating string within field theory [1] and other problems jl5] . Locally speaking, 
these problems correspond with Lagrangian and Hamiltonian functions that do not depend on the 
base coordinates, usually denoted by [t^ , . . . ,t^). When a dependence on the base coordinates 
exists, the /c-cosymplectic formalism is necessary. In other words the A;-cosymplectic formalism 
is the generalization of the cosymplectic formalism used in non- autonomous mechanics to field 
theories [T^fTB]. 

However, the control equation ([1]) under study is of second order. Hence we need to extend 
the A:-symplectic formalism for optimal control problems developed in this paper to implicit control 
differential equations. In this framework we transform the equation ([T|) into a first order one in such 
a way that we obtain an implicit equation. 

The paper is organized as follows: In Section [2] we define the setting to describe optimal control 
problems governed by an explicit first-order partial differential equation using the A:-symplectic 
formalism. A parallel between this formulation and the geometric description of optimal control 
problems governed by ordinary differential equations is considered in order to stress the similarities 
and differences between both problems. 

One of the key points to prove Pontryagin's Maximum Principle consists of extending in a 
suitable way the control system by adding new coordinates which contain the information related 
to the cost function, the so called extended system. In this paper, the /c-symplectic formalism for 
the extended optimal control problems only works under some particular assumptions on the cost 
function, which turn out to include the most typical cost functions in the literature. To prove 
Pontryagin's Maximum Principle on fe-symplectic formalism in Section [2. II we define the elementary 
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perturbation vectors on that formalism. 

After this first approach to tackle optimal control problems, we are going to consider in Section [3] 
the unified Skinner- Rusk formalism for fc-cosymplectic implicit dynamical systems. Following the 
lines of [2], we adapt the above formalism to describe a novel unified formalism for optimal control 
problems governed by an implicit partial differential equation in Section SI This generalized unified 
formalism will allow us to consider interesting problems associated with higher order control partial 
differential equations, in particular the problem that has motivated our study, see Section [5j In this 
last section, we consider the control partial differential equation that models the orientation of a 
bipolar molecule in the plane studied in [6], as described at the beginning, when a control-quadratic 
cost function is considered. Hence, the use of A;-symplectic formalism and all the generalizations 
we have made in the previous sections to deal with optimal control problems on partial differential 
equations are fully justified. 

In the sequel, unless otherwise stated, all the manifolds are real, second countable and C°^. The 
maps are assumed to be also C°°. Sum over all repeated indices is understood. 



2 /c-symplectic formalism for optimal control problems governed 
by an explicit first-order partial diff^erential equation 

We first recall briefly the essential definitions and notations in the /c-symplectic formalism. Let 
Q be a n-dimensional manifold and tq : TQ — >■ Q be the natural tangent bundle projection. The 
k-tangent bundle or the bundle of -velocities of Q, denoted by T^Q, is the Whitney sum of k 
copies of the tangent bundle TQ, that is, 

TlQ = TQ® ®TQ. (2) 
The elements of T^Q are /c-tuple (vi^, . . . , Vkg) of tangent vectors on Q at the same point q G Q. 
The canonical projection Tq: T^Q ^ Q is defined as follows 

rQ{vi^,---,Vk^) =q. 

If (y, {q^)) is a local chart on Q, then it induces a local chart {T^V, on T^Q, where T^V = 

A k-vector field on Q is a section X: Q ^ canonical projection Tq. Hence a k- 

vector field X defines a family of k ordinary vector fields {Xi, . . . , X).} on Q through the canonical 
projections Tq : T^Q — > TQ onto the A-th component of T^Q, that is. 



T0^ivU,---,Vk„) =VA 
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where A = 1, . . . ,k. Note that Xa = Tq^ o X. 



Q 

An integral section o/ X is a map o" : M'^ Q, t = (t-^, . . . , t'^) — > a{t) such that 



rpi I da da\ 



We introduce now the notion of control system in the /c-symplectic formalism. Consider a control 
set U We need the notion of a /c-vector field X defined along the projection vri : Q x {/ — >• Q. 

Such a /c-vector field is defined by making the following diagram commutative: 

J C M'^ QxU ■ 
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where J is a subset of . 



An integral section of a k-vector field X defined along the projection tti : Q x U — )■ Q is a map 
= {a,u): J CR'' ^ Q X U such that 



(Ti)Ko,^) = Xo, 



or in other terms, 



d 



V,...,j.)(7ri 0^)^= XAm\ . . . , t^)) = XA{a{t\ t''),u{t\ t")), 



dt^ 

Xa ° 4> for A = 1, . . . ,k. 



equivalently (tti c 

Let F: Q X U ^Rbea regular enough map. Such a function, which is usually called the cost 
function in the literature, allows us to define the functional 



n 



Dom I 
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where dS = dt^ A • • • A dt'^, i.e. the usual volume form in M.^. 

From now on, we assume that Dom (p = Ii x ■ ■ ■ x If, = [tQ,t^] x ■ ■ ■ x [t^, t'j] = : I. 

Before stating the optimal control problem on A;-symplectic formalism, we remind here the clas- 
sical optimal control problem with cost function G : Q x U M. 

Statement 1 (Optimal control problem, OCP). Given {Q, U, X, G, I). Find a curve (7, u): I C 
R ^ Q X U joining the points xq and xj va. Q such that 

(i) it is an integral curve of the vector field X defined along the projection tti: Q x U ^ Q, i.e. 

j{t) = x{j{t),u{t)y, 

(ii) it minimizes the functional Jj G{'y{t),u{t))dt among all the integral curves (7, tt) of X on QxU 
joining xq and Xf. 

Statement 2 (/c-symplectic optimal control problem, /c-OCP). Given (Q, C^, X, F, I). Find a 
map (f) = {a,u) : I = Ii X ■ ■ ■ X Ik C M.'^ ^ Q X U passing through the points qo and Qf in Q such 
that 

(i) it is an integral section of the /c- vector field X = {Xi, . . . ,Xk) defined along the projection 
TTi: Q X U ^ Q, i.e. 



TiKo0) = Xo0, i.e. ■^it) = X\{m) = X\Kt),uit)), 



(4) 



where t = {t^,...,t^) G h 



X 4; 



(ii) it minimizes the functional Jj^^ F{(j){t))d^t among all the integral sections ^ of X on 
Q X U passing through qo and qj, where d'^t = dt^ A ■ ■ ■ A dt'^ . 

Let us compare the frameworks in the traditional optimal control problems and in the k- 
symplectic optimal control problems. 



Classical OCP for ODE 
TQ 



fc-symplectic OCP 



QxU 



X 



Q 



(7,w 
I C 
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QxU 

(j) = (cr, u) 

1 = h X ■ ■ ■ X Ik C 
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Observe that the classical OCP has associated a problem of explicit ordinary differential equa- 
tions, whereas the equations in the A:-symplectic optimal control problem are explicit partial differ- 
ential equations. 

Since late fifties the most efficient tool to solve optimal control problem is Pontryagin's Maximum 
Principle, which provides us with necessary conditions for optimality [16]. One of the key points to 
prove that Principle for classical optimal control theory consists of extending the control system in 
a suitable way. To be more precice, Q is extended to the manifold Q = R x Q with local coordinates 
X = and the corresponding extended vector field is given by 

- ^ d 

X(x,u) = G(x,u)^—^ +X(x,u). 

OX^ {x,u) 

Note that the system of ordinary differential equations which determines the integral curves of X 

dx^ 

can be decoupled in the following sense: we first integrate = X^{x,u) and then we have 

x°{t)= f G{j{s),u{s))ds, 
Jto 

for any t G / = [to,ii]. 

Unfortunately, in order to extend coherently the optimal control problem on /c-symplectic for- 
malism we need some extra assumptions on the cost function. 

Assumption 1. The Lie derivative of the cost function with respect to each Xa is zero, that 
is, Lx^F = 0. 

Assumption 2. The control functions n: I — )• C/ are locally constants. 

We can justify these assumptions as follows. In a first try to extend the control system we will 
add k new variables {q^^ , • • • , (p^) such that for every A= 1, . . . , /e 

- F. (5) 



dqOA 



0, for B^A. (6) 



Once we have an integral section of X, we integrate Q and obtain 

q^^[t\...,t^,...,t'']{t^)= / F{q\...,q^,u\...,u'){t\...,s,...,t'^)ds. 



The equation ([6|) is satisfied by q^^ for B ^ A'li 

dqOA ^ rt^ /QF dq' dF du''\ ^ ^ 

Note that if both assumptions are satisfied the equations will be immediately satisfied. Thus these 
assumptions are necessary to guarantee the compatibility of the system of partial differential equa- 
tions when we extend the control system. 

Moreover, having in mind [U] Assumption 2 is reasonable when dealing with control systems. 
These two assumptions include the most typical cost functions considered in optimal control prob- 
lems such as control-quadratic, constant function 1 (that is, time optimal), etc. Hence these as- 
sumptions do not impose great restrictions according to the literature. 

Under the above assumptions, let us consider now the extended /c-symplectic optimal control 
problem. In order to preserve the same philosophy as in classical control theory, we will have to add 
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k new coordinates {q^^){B=i,...,k}- Then the extended manifold in fc-symplectic formahsm is given 
hy Q = M.^ X Q. If X = {Xi, . . . , Xk) is the A:-vector field on Q, then the extended k-vector field X 
on Q is given by {Xi, . . . , X^) where 

^^ = ^^^5^ + ^^ = ^9^ + ^^' A = l,...,k, (7) 

where is the Kronecker's delta and F is the cost function. 

As in classical optimal control theory, we can first integrate = X\{q,u) and obtain (j) = 
{a,u):I^QxU. Then, 



gOA(t) = / F(a(t),n(t))dt^, for every A = l,. 



,k. 



Because of ([6]) q^^ is constant when we fix t^. On the other hand, due to Assumption 1 andOq*^-* 
is constant along integral curves of Xa_ for every A = 1, . . . ,k. 

Statement 3 (Extended A;-symplectic optimal control problem). Given (Q, [/, X, F, I). Find 
a map (f) = {a,u) : I C M.^ ^ Q x U passing through the points (0, qo) in Q and qf in Q such that 

(i) it is an integral section of the /c- vector field X = {Xi, . . . ,XjS) defined along the projection 
vfi : Q X [/ — )■ Q, i.e. locally 

Ti{^,o4>) = Xo4>, i.e. -^(t) = F{4>{tM, ^{t)=X\{m), 

where t = {t^ , . . . ,t^) £ Ii x ■ ■ ■ x 1^, 6^ is the Kronecker's delta, for every A,B = 1, . . . ,k; 
i = l,...,n; 



(ii) it minimizes each functional 



t)= / Fm\...^f,...,t^))ds, (8) 



for A = 1, . . . ,k, among all the integrals sections </) of X on Q x [/ passing through qo and qf 
such that (j) = ttqxU ° 4> foi^ T^Qxii ■ Q x U ^ Q x U. 

Remark 4. Note that if the projection of an integral section (p : Ii x ■ ■ ■ x 1^ —?■ Q x U of X to on 
Q X U minimizes each functional in ([8|), then the projection of the integral section (j) to (f) on Q x U 
minimizes the functional 



m= / (Fo,^)d^t (9) 

since the order of integration does not matter. 

Hence, in contrast with classical optimal control theory, in /c-symplectic formalism the extended 
optimal control problem and the optimal control problem are not equivalent. However, solutions to 
the extended problem in Statement [3] are also solutions to the original fc-symplectic optimal control 
problem in Statement [2l As we will see later on, the adapted version of Pontryagin's Maximum 
Principle in fc-symplectic formalism provides us with necessary conditions for optimality of the k 
functionals in ([8]) for those cost functions satisfying assumptions 1 and 2. 

As mentioned above the trajectories that minimize ^ also minimize ([9]), but not necessarily in 
the other way around. Remember that to minimize a multiple integral does not imply that every 
simple integral involved is minimized. Thus, the necessary conditions for optimality described in 
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Section [2. II in the A:-symplectic version of Pontryagin's Maximum Principle are more restrictive than 
the traditional necessary conditions for optimality in [3l [16] . 

The elements of extended optimal control problems in classical formalism and A;-symplectic 
formalism are summarized in the following diagrams: 



Classical extended OCP for ODE 



fc-symplectic extended OCP 





• • X 4 c 

In order to state Pontryagin's Maximum Principle we need a Hamiltonian problem associated 
with each of the extended optimal control problems. It is important to remark here that this 
Hamiltonian problem is not equivalent to the optimal control problems in the classical formalism 
neither in the /c-symplectic formalism. 

The Hamiltonian for the extended optimal control problem in classical theory is given by 
H:T*QxU^R, 

H{p,u) = {p,X{x,u)) =poG{x,u) +piX'{x,u). 

For each control u, the Hamiltonian vector field = Xh{-,u) satisfies the following Hamilton's 
equation 



di7^"> 



where oj is the canonical symplectic structure on T*Q. Locally u = dx^ A dpo + dx* A dpi in natural 
local coordinates {x^ ,x^ ,pQ,pi) in T*Q. 

For the extended A;-symplectic optimal control problem we consider k Hamiltonian functions 
Ha - {TlYQxU ^R defined as follows 

fc n 
//a(p, u) = {p^,XA{q, u)) = Vt^iq, n)6^ + Y,pfx{{q, u) 



B=l 



p^^F{q,u)+Y,pfx{iQ,u), 



(10) 



in natural local coordinates (g°S . . . ,q^'' ,q^ , . . . ,q"' , {pq^, . . . ,pQ^,pf, . . . ,Pn)A=i,...,k) for {T^)*Q. 



For each control u, the Hamiltonian A;-vector field X' 
following family of equations 



Xj**"* , • • • , ^fc*"* ) must satisfy the 



^o,M^^ = dH^^ for every A = 1, . . . , k. 



(11) 



The canonical /c-symplectic structure on {T^)*Q is given by (cji, . . . , w^) where ua = (7r^)*w, 
TT^ : [T^YQ T*Q is the projection onto the Ath-copy and uj is the canonical symplectic structure 
on r*Q. Locally loa = dq°^ A dp^^ + dq' A dpf. 

If for each control the Hamiltonian /c-vector field X* is solution to ()lip . then it is solution to 
the following Hamilton-De Bonder- Weyl equations 



A=l 



{u} 
A 



dH, 



(12) 



A=l 



K.A=1 
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associated with the Hamiltonian H: {T^)*Q x [/ — > M given by 

k k 
H(p,n) = J2Ha{p,^) = Y,{p^,XA{q,u)), 
A=l A=l 

where p G (^a:)|^'3- By the superposition principle all the solutions of (fTTj) are solutions to (fT2]) 
because both systems are linear in the momenta. However, ()12p has more solutions apart from the 
ones coming from (jlip . In fact, for every A € {1, . . . ,k} the Ath vector field of the Hamiltonian 

fc- vector field X* = (x^, . . . jX^j is locally expressed as follows 



From (jlip we obtain 



A _ A dF . dX\ (13) 



dq^ dq^ 



A=i A=i 



for every A = 1, . . . ,k. Note that the Hamiltonian A:-vector field X* = yX^, . . . , X^^j is not com- 
pletely determined because the following functions 

{Ya)1, {YA)f (14) 

remain undetermined for C A and for every A = 1, . . . ,k. 

On the other hand, from (jl2p the Hamiltonian /c-vector field X* = (^Xl, . . . , X^^ must satisfy 

{Ya)'"> = X"/ = FSl {Ya)H,=0. {Ya)' = X\, 

By comparing (jl3p and (|15p it is clear that all the solutions to (|lip are also solution to (jl2p . 
but not in the other way around. Neither the Hamiltonian A:-vector field X* " solution to (jlip . nor 
the Hamiltonian /c-vector field X* solution to the Hamilton-De Bonder- Weyl equations are fully 
determined. For the first one, the functions in p4p remain undetermined. For the second one, the 
functions in (|14p remain undetermined and maybe some of the {Ya)^ involved in ()15p . 

However, we can reduce in an intrinsic way the number of functions that remain undetermined in 
the above mentioned Hamiltonian /c-vector fields X**"^ in such a way that the Hamiltonian /c-vector 
field X*^"^ solution to ^ is fully determined. Note that {T^)*Q = {T^T{M.^ x Q) ~ {T^yR^ x 
{T^)*Q, which has two natural projections pr^^ and pr2 from {T^)*Q to {T^yU.^ and {T^)*Q, respec- 
tively. Consider now the canonical projections tTq'*^: {TI)*Q T*Q and vrjf : {T^)*R'' T*R'' to 
the Cth component of {T^)*Q and {T^)*R^, respectively. 

The conditions 

T(<f opri) (Xl) = 0, (16) 

T(^n'fopr,)(^XX) = 0, (17) 

for every C ^ A imply locally that {Ya)qb = and (^)^ = for C A and for every C A. 

Under conditions ()16p . ()17p . given an initial condition /3o in {T^)*Q there exists a unique integral 
section /3 : /i x • • • x — > {T^)*Q of the Hamiltonian /c- vector field solution to (jlip . It is clear from 
the local equations (|13p that once A is fixed, p^ does not appear in the set of equations in (jlip 
associated with A and only p^^s appear. 
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2.1 Elementary perturbation vectors and Pontryagin's Maximum Principle on 
/c-symplectic formalism 

Now let us introduce the notion of elementary perturbation in fc-symplectic formalism that allows 
us to define later the /c-symplectic tangent perturbation cones. These elements are essential to prove 
the A;-symplectic Pontryagin's Maximum Principle, Theorem [8l 

First fix a surface (a, n) : /i x • • • x — t- Q x [/. Let vr^ be a 3-tuple {va, Ia, ua} where va, 
G M and ua ^ U C M} . The Ath-elementary perturbation of the control u is defined as 
follows ^ 

u[7r%]{t\ = { ... ^ ^fc)^ elsewhere. ^^^^ 

Associated to this control ti[7r^], the mapping ?[7r^] : Ii x • • • x — > Q is the integral section of 
the /c-vector field X.^'^^^a]} with initial condition (tj, . . . ,tQ,a{tl, . . . ,1^)). 

Given e > 0, define the map 

: /i X • • • X 4 X [0, e] Q 

{t,s) I — > f^^{t,s) = a[7r%]{t). 

For every t G /i x • • • x 4, tp^^ : [0, e] ^ Q is given by <^*^(s) = ¥5^^(t, s). The curve ip\^ depends 
continuously on s and on it a = {rA,lA,UA}- 

Prom 5[7r^] we can define a curve as follows 
a[7r%]it\...,t^,...,t>'): Ia Q 

This curve is an integral curve of A^"''^'*^^ with initial condition {tQ,a{tQ, . . . ,to))- 

Proposition 5. Let va G Ia- If '^[''^a\ elementary perturbation of u specified by the data 

T^A = {'''Aj^AjUa} , then the curve (p\^ is differentiahle at s = and its tangent vector is 

[XA{a{t\...,rA,...,t''),UA)-XA{d{t\...,rA,...,t''),u{t\...,rA,...,t''))]lA =: v[7^a] (19) 
for fixed . . . , r^, . . . , i'') G /i x • • • x 4. 

Proof. In local coordinates (q^"^ , ■ ■ ■ , <f'' ,Q^,---,q"')^otQ, note that 

q' o d[7r%]{t\ ...,tA,.. .,t^){rA) - q' o dWA]{t\ . . . . . . ,t^){t^) 

X\{a[iT\\{t\ ...,tA,.. .,t>'){t),u[7r\]{t\ ...,t,... t'')) dt 



for every i G {Oi, . . . ,Ofc, 1, . . . n). 

To compute the derivative of 93*^ at s = with t = (t^, . . . , rA, ■ ■ ■ ,t^) we use the definition of 
the derivative: 



_d 
d^ 



W °¥',r^)('S) = lim' 



lim ° ...,rA,...,t^)-q'o a{t\ . . . , r^, . . . , t^) 

s^O S 



lim 



7;/ X\{d[i,\]{t\ ...,tA,.. .,t'){t),u[7r%]it\ ...,t,... t>^)) dt 
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/;/ X),{a{t\ ...,t^,.. .,t''){t),n{t\. ..,t,. ..t''))di 



lim 

s-i>0 



/;;_. ^ x\{a{t\ m), u{t\ ...,t,... t')) dt 



c. 



Let us use now the following equation 

/ X{j{h),u{h))dh = sX{j{t),u{t)) + o{s), (20) 

Jt-s 

in the above formula having in mind that o(s) tends to when s tends to 0. Then, 
' X'^{d[7r%]{t\ ...,t^,.. .,t^){rA),UA)lAS 



C = lim 



X\{a{t\ ...,tA,..., t>'){rA),u{t\. ..,rA,... t^))lAS + o{s) ' 



s 



= hm (x\{d['K\]{t\ .,t''){rA),UA)lA 
s— >-0 \ 

-X\{d{t\ ...,t^,..., t'')irA),u{t\ ...,rA,... t''))lA) 

= [Xi(a(t\...,rA,...,i'=),^XA)-X^(?(t\...,rA,...,t^),n(t\...,rA,...,t'=))]/A := 
for each i G {Oi, Ofc, 1, ... n}. ■ 

Note that the tangent vector in Proposition 12.11 is in ( r'^'^iTlQ)] = T^/.i ik\Q. The 

\ Q /5(ti,...,tg) '^(hi'-'h-,)^ 

vector vItta] is called the Ath-elementary perturbation vector associated to the perturba- 
tion data TTA = {rA,lA,UA}- It is also called an Ath-perturbation vector of class I. 

Following the same lines as in [3] we can define the associated Ath-perturbation vector obtained 
from c different Ath-perturbation data ttai, ■ ■ ■ ,t^Ac with different and/or same perturbation time 

At each copy of the tangent bundle in the fc-tangent bundle T^Q, we construct an Ath-tangent 
perturbation cone 

, a<r<t 



= conv U (21) 

\a<T<t ^ J * J 

where denotes the set of Ath-elementary perturbation vectors at r, ( ^^^^^ 1 is the push- 

forward of the flow of A^ with ^(r) as initial condition at time r, convM^ denotes the closure 
of the convex hull of the set W . 

Remark 6. If the controls are only measurable and bounded, as usually assumed in control theory, 
all the perturbations and geometric elements such as vectors, cones, etc. that appear in the paper 
are only defined at Lebesgue times where the equality (pJj) is satisfied. 

The definition of k different perturbation cones in ()2ip implies that the perturbation data as- 
sociated with different Ath. copies are not mixed. As proved in [3l Proposition 3.12], the following 
result is true for the cones for every A = 1^ . . . ,k. 
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Proposition 7. Let G [to^tf]- If v is a nonzero vector in the interior of K^, then there exists 
e > such that for every s € (0, e) there are s' > and a perturbation of the control u[ita] such that 

a[7r%]{t\ . . . , t"^, . . . , = a{t\ . . . , t'^, . . . , t^)(t^) + s'v. 



This proposition is essential to prove Pontryagin's Maximum Principle in /c-symplectic formalism. 

Theorem 8 ( A;-symplectic Pontryagin's Maximum Principle) . If (p* = (a* , u* ) : /i x • • • x — )• Q x [/ 
is a solution of the extended k-symplectic optimal control problem {Q,U,Ji.,F,I), Statement\^ such 
that F satisfies assumptions 1 and 2, then there exists (/?, u) : /i x • • • x 1^ — )• {T^)*Q x U along a* 
such that 

(i) (vr"^ o /3, u) along a* is a solution of (jlip for each A = \, . . . ,k; 

(a) the Hamiltonian Ha'- {T^)*Q x [/ M in (jlOp along the optimal integral section is equal to 
the supremum of Ha over the controls almost everywhere; 

(Hi) the supremum of the Hamiltonian Ha'- {T^)*Q x U — t- M in (jlOp along the optimal integral 
section is constant almost everywhere; 

(iv) /5^(t) 7^ G T^,(t)(5 for each t e h X - - - X Ik and for every A = 1, . . . ,k; 

(v) /S,^ (t), . . . , /3,^ (t) are constant and (3^^ is non-positive for every A = 1, . . . ,k . 

Proof. As (a*,u*) is a solution of the extended fc-symplectic optimal control problem, if r G /i x 
• • • X Jfc, for every initial condition /?,- in (T^)*Q there exists a unique curve /3 in {T^)*Q satisfying 
the k equations in (|lip and the initial condition. 

As in the classical Pontryagin's Maximum Principle the initial condition must be conveniently 
chosen so that the rest of conditions in the theorem are fulfilled. 

For each ^4 G {1, . . . , k}, consider the Ath-tangent perturbation cone KA C T-^*u,)Q and the 

A ^ 
vector (0, . . . , —1, . . . , 0, 0, .".,0) in T-^*(tf)Q that indicates the decreasing direction of the coordinate 

q'^it) = JtA Fia*it\---,h,...,t''),u*it\...,h,...,t>'))dh. 

A 

Observe that if (0, . . . , — 1, . . . , 0, 0, 0) was in the interior of Kix, then there would exist an 

^ / 

Ath-perturbation data tta = {rA, Ia,ua} such that (a[7rA], ^^[vrA]) passes through the same points on 

Q as a* = TT2oa* ,hut q^'^lal'KAjlitf) < q^'^[d*]{tf). This is a contradiction with the fact that (5*, ti*) 

A 

is a solution of the extended k-symplectic optimal control problem. Hence (0, . . . , — 1, . . . , 0, 0, 0) 
cannot be in the interior of K^. 

Thus, there exists G T^*(^tf)Q such that 

(^;i,(0,...,-l,...,0,0,.'?.,0)) > 0, (22) 
0fA,v[TTA]) < V%A]ei^tj- (23) 

Condition ()22p implies that (3^^ < 0. Let us explicitly write for each A G {1, . . . ,k} the equations 
for the integral curves of X\ that satisfy equation (fTTj) : 



F61. -1^ = 0, 



dt^ " ^' dt^ 

^ = ^ = _ A OF ,dX\ 

Q^A A^ Q^A Qqi I'OAQqi l^j Q^i ' 
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for B = 1, . . . ,k and i,j = l,...,n. Note that there are no equations for pf with B A. Hence the 
momenta whose coordinates are pf remain undetermined. They will be determined by solving (jlip 
with A = B. Given an initial condition in the Ath-copy of T^^^^^^Q, we just solve the equations in 
the fiber for p^. 

If = G T^m^)Q, then the solution to (HU in the fiber will be zero along a* because of the 

linearity of the differential equation in the momenta. If the momenta is zero, it does not provide 
us with any information related to the separation condition. Hence, P^{t) 7^ G ^^•(t)0 for every 
t € I. 

As in the classical Pontryagin's Maximum Principle, condition ()23p and the definition of the Ath- 
elementary perturbation vector in p9|) prove the condition about the supremum of the Hamiltonian 
Ha- {T^)*Q xU ^ M in (|10p over the controls for each A G {!,..., A;}. The constancy of the 
supremum of the Hamiltonian over the controls is proved analytically, analogously to the classical 
Pontryagin's Maximum Principle, see [3l \W\ [T9] for more details. 

Prom equations ([TT]) we deduce that for each A € {1, . . . ,k}, (3q^ is constant for every B G 
{1, . . . ,k} along the optimal integral section (p* = (a* ,u*). ■ 



3 Application of unified formalism for /c-cosymplectic to implicit 
PDEs 

We are going to apply the unified Skinner- Rusk formalism for fc-cosymplectic field theories developed 
in [n\ Section 4] to the dynamics description for systems given by implicit partial differential 
equations. This will be very useful to develop Sections H] and [5] so that physical examples associated 
with higher order control partial differential equations fit in the approach considered in this paper. 

The Whitney sum T^Q (T^)* Q has natural bundle structures over T^Q and (T^Y Q. The 
suitable bundle to describe non-autonomous dynamical systems governed by partial differential 
equations is W: = M'^ x [T^Q ® [T^)* Q). Local coordinates for W are {t^ , , v\, pf) . Let us 
denote by pr^ : M'^ x {T^Q © (T^)* Q) ^ ]R^ pra : x {T^Q © (T^)* Q) T^Q and prg : M'^ x 
(T^Q ® (T^Y Q) ~^ i'^k)* Q local projections into the first, second and third factor of W, 
respectively. Locally, 

pT,{t^,q\v\,pf) = t^, 
pr2it^,q\v\,pf) = iq\v\) = iq,v), 
V^S^ ,q\v\,pt) = {q\pt) = {q,p). 

Let (dt^, . . . , dt^) and {ui, . . . ,ujk) be the canonical forms on M.^ x (T^)*Q. We denote by 
. . . , i?^ and (f^i, . . . , Ofc) the pullback by pr^ and prg of these forms to M'^ x {T^Q © (T^)* Q) , 
that is, i?^ = (pri)*(dt^) and Qa = {pt3)*{uja) for 1 < A < A;. Locally, 

= dt^, Qa = dq' A dpf. (24) 

The coupling function C on R'' X [T^Q © [TlY Q) is defined as follows 

C: R'' X tIQ © {Tiy Q ^ R 

k k 

(t,V„p,) ^ J^P^(^A) = J^(^'^^^a)- 

A=l A=l 

Given a Lagrangian function L on M'^ x T^Q, the Hamiltonian function H on R^ x (T^Q © {T^)* Q) 
is defined as follows 

H = C- (pri X pr2)*L. (25) 
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Locally, H{t, q\v\,pf) = pfv\ - L(t, q\ v\). 

The problem in the Skinner-Rusk formalism for /c-cosymplectic field theories consists of finding 
integral sections (p: R'' x {T^Q ® {TD* Q) of an integrable /c-vector field Z = {Zi, ... , Z^) on 

R^ X {TlQ e {tIYQ) such that: 

jZ-^z,^A = dH-j2^^\ ■^zJ^ = Sl (26) 
A=l A=l 

See [T7] for more details. 

After summarizing briefly the Skinner-Rusk formalism for /c-cosymplectic field theories, here we 
are interested in adapting it to find the dynamics of systems described by implicit partial differential 
equations. An implicit dynamical system (L, M) is described by the submanifold 

M = {{t'',q\v\) € M^^ X {Tl)Q | ^^{t'' ,q\v\) = 0, 1 < q < s}, 

of R^ X TlQ, where d^^ A • • • A d^''* / 0, and a Lagrangian function L € C°°{M). This submanifold 
M of M'' X TlQ can be naturally embedded by l^' : M ^R^ x TIQ. 

In order to adapt the above formalism to this kind of dynamical systems, we must define the 
fc-symplectic implicit bundle W*^ = M xq {TI)*Q and the corresponding canonical immersion 

iM . y^M ^y^. ^^k ^ (^IQ ^ ^^1^ * ^ (27) 

Now we can consider the pullback of the coupling function and the canonical forms on W to W^'^ : 

c>^" = (i^^)*(C), = (i^^)*(^^), nT = i^''n^A). 

Let : W*^ ^ M be the natural projection, we define the Hamiltonian function Hy^M : 
R as follows 

H^n = C^" - (pf )*L. 

Analogously to (j26p . the problem of describing the dynamics of (L, A/) consists of finding the 
integral sections M'^ — )• W^^ of an integrable A;- vector field Z = {Zi, . . . , Z^) on W*^ such that 



j^i^.^^r =d//w--E^<-' i^.dt^=5f- (28) 

A=l A=l 



Or equivalently, the problems consists of finding the integral sections M*^ — )• W of an integrable 
fc-vector field Z = (Zi, . . . , Z^) on W = M^' x {tIQ [TI]* Q) such that 

Ei^.f^^ = d^wM-i:^fr^^ + ^"dvI."-A„^^^^ iz.^^ = 6l (29) 
A=l ^=1 A=l 



This equation is obtained from ()28p by rewritting the equations on W so that the constraints = 
in (j27p must be added to the equation in a suitable way. 

The Ath vector field Z^ on W is locally given by 
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Prom ([29]) we first have {Za)^ = 1, (^a)? = fov B A. Moreover, 
A=l A=l A=l 

= {ZATdpt - - v\dpt - pfdv\ + l^dg^ 



Thus, 



[ZaT = v\, (30) 



By also imposing the conditions (jl6p and ()17p in the /c-vector field on W, we have 
with 



^ _ ^ _ 9^ -A , sA _ ^ _ X ^ 

If Z is a solution of (j29p . then we must start a constraint algorithm in the sense of [?]• To 
be more precise, each Za must be tangent to the submanifold Ml contained in W^^ and defined 
by (|32p . That is, the following tangency conditions must be satisfied 

= ZAi-^n, (33) 
= pf- — + , (34) 

V OV^ OVg J 

on Ml for every A = 1, . . . ,k. Depending on the particular examples, some components will be 
determined and the constraint algorithm must proceed until stabilization. 

Remark 9. For non-autonomous fc-symplectic explicit dynamical systems the manifold M is defined 
by constraints locally given by \I'°(t, g% f^) = v^a ~ ~ ^- Hence the above process can be 

used for this kind of dynamical systems. 



4 Unified formalism for optimal control problems governed by an 
implicit partial diff^erential equation 

We extend now the unified formalism for implicit control systems developed in [21 Section 4] to 
optimal control problems whose dynamics is given by implicit control partial differential equations, 
instead of just explicit control partial differential equations as developed in Section [2j The problem 
consists of finding the solutions to optimal control problems governed by an implicit partial differen- 
tial equation by taking advantage of the unified formalism developed in this section. See Section [S] 
for a particular problem where this unified formalism is used. 
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Let C be the control bundle with natural coordinates (t"^, g*, ti*^). In contrast to the explicit 
description of control partial differential equations in Section [21 let us consider now the case where 
the control partial differential equations are given implicitly by the following submanifold 

Mc = {(^^,n^g^^;^) e C XqT^Q \ g^ ^;^) = 0, 1 < a < s} 

of CxqT^Q, where d'l'^ A • • ■ A d^''' ^ 0. There exists a natural embedding l^c Mc ^ C xqT^Q- 
Then the implicit optimal control problem under consideration is determined by (L,Mc), where 
L G C°°{Mc) is a Lagrangian function. 

Let us define now the /c-symplectic implicit control bundle W*'^"^ = Mc Xg [T^)*Q which is a 
submanifold of C x^ky^gW = C x^ky^q (M*"' x {T^Q © {T^)*Q)). Then we have, respectively, the 
canonical immersion and the natural projection: 

iMc.y^Mc ^Cx^.^qW, an^-.Cx^k^gW^W. 

Now we can consider the pullback of the coupling function in Section [3] and the canonical forms on 
W to W*^c: 

C>^«^ = (^^ o i'^cy(^c), = (^w o i^^)*(0^), = (aw o i^^)*(^^) 

Let pf^^ : yV^^^ Mc be the natural projection, the Hamiltonian function Hy^M^y : W*^"^ — )• M 
is defined as follows 

The dynamics of the optimal control problem (L, Mc) is determined by the solutions of the equations 

Eiz..(f^r")=0, iz.KMa=Sl (35) 
A=l 

for a A;- vector field Z = (Zi, . . . , Zfe) on W^^c^. 

In order to work in local coordinates we need the following proposition whose proof is straight- 
forward. 

Proposition 10. For a given w G W^^c-^ ^/jg following conditions are equivalent: 

(1) There exists a k -vector field E (r^)«;W*^^ verifying that 

k 

E 0^"^((Z^)^, {YaU = , for every Y,^ G {T^^W^'^ . 
A=l 

(2) There exists a k -vector field S (T^)w{C x^ky.Q W) verifying that 

(i) G {tD^w'''^, 

(^^) EA=ii(Z.)J^w(^^A))«,G((Ti).W^-^^)0 . 

As a consequence of this last proposition, we can obtain the implicit optimal control equations 
using condition (2) in Proposition [10] as follows: there exists a fe-vector field Z on C Xjgfc^g W such 
that 

(i) Z is tangent to W'^'^ ; 

(ii) the 1-form i(ZA)('7vv(^A)) is null on the /c-vector fields tangent to W^^^ . 
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As W^c- = McXq {TI)*Q and the constraints are = 0; then there exist G C°°(C XjjfexQ W), 
to be determined, such that 

J2Hz^){a*,^{^A))\wM, = d/7^Me -^^^^ + A„d^"-A„^^9^ \^Ma. (36) 
A=l ^ ^ 

As usual, the undetermined functions A^'s are called Lagrange multipliers. 

Now using coordinates {t^ , u°- , q\ v\, pf) in C x^k ^g W, we look for k vector fields 

= (Z.)f ^ + + + + (F.)f ^ . 

where (Za)^ , (-^a)", (-^a)*) (-^a)^) (^A)f are unknown functions on W^'^'^ verifying the equation 



k / k \ 

= J^iz^(dg^Adpf)-d J^(pM)-L(t,u,<?,v) + 
A=l \A=1 / 



-A„d*° + A„^dt^ 

Note that from ([35]) we have {Za)^ = 1 and {Za)^ = Oior Ay^ B. Moreover, 

(^a)' = (38) 
Pf = WT-^o.jrT^^ (39) 

together with the tangency conditions 



= ^A(^")lw-c = { iZA)t^ + {ZaY^ + (ZaT^ + (ZaTb^ ) (41) 



for every A = 1, . . . , k. Imposing the conditions ((T6]), ()17p . we know that (Ia)^ = for A 7^ i?. From 
here we can start a constraint algorithm in the sense of ^ as follows: the tangency conditions with 
respect to the constraints ([39|) and (P0|) obtained from equation ([36|) give the following equations on 
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y^Mc. 



9L , dib"^ 



dv 



B 



dv 



B 



dt^dv 



+ K 



dqidv 



B 



{Za)^ + \ 



du°-dv 



B 



du°-dv^^ 



[ZaT: 



Za 



+ 



dt^du"- du^dw 



(Za] 



B 



Aq, 



dq^dw 



(ZaYb - Xc 



dvPdi 



If the square matrix of size k + nk 

( 



Aq, 



+ A, 



a2$« 

92 

dvPdvQ 



— Aq, 

+ A, 



a2^a \ 

dv^gdu"- 



has maximum rank, (.Z^^)^ and (ZaYb determined in terms of (Ya)^^, which must satisfy the 
condition pzp coming from (j36p . The algorithm continues until stabilization. 



5 Example: Orientation of a bipolar molecule in the plane by 
means of two external fields 

Let us consider now the control partial differential equation studied in [6l Section 8]: 
.d^it,9) ( 



-y +ui{t)cose^>{t,e) + u2{t)sine^{t,e)y (42) 

where ^ is an element in a Hilbert space taking values on the complex and ui, U2 take values in M. 

Let us rewrite the problem according to Section [H This equation fits in 2-symplectic formalism 
where = t and = 9. Note that ()42p is a partial differential equation on the complex numbers. 
Hence let us rename 

q^ = Re^*, q^ = Im^. 



In order to rewrite ()42p as an implicit partial differential equation we work on a 6-dimensional 
manifold Q with local coordinates 

12 3 dq^ i dq"^ 5 6 9q'^ 

'^'^'^ =^'^ =^''^ =^ 

to transform the second order partial differential equation into first order partial differential equa- 
tions. 

The local coordinates for C xg TgQ are {t^,t'^,u^,u'^,q^,v\,V2). Note that apart from (j42p we 
also know that 



Vl 



' dq^ dq^ d'^q^ d'^q^ d'^q^ d'^q^\ f dq^ dq^ d'^q^ d'^c^ d'^q^ d'^q'^\ 
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Hence (03]) determines some relationships between some coordinates of vi and V2- Equations ()42p 
and ()43p . determine a submanifold Mc of C XgTgQ implicitly defined by the following constraints: 



^1 






xI/4 


— ^2 




= —q'^ — ^2 + uiq^ cos 9 + n2(7^ sin 


^2 








= 




= 9^ — t'l + ^1*?^ cos + ti2g^ sin 6 


^3 






^6 


— ^2 


-vf, 





A general 2- vector field Z on C Xi[j2xq (K^ x (Tg^Q © {T2TQ)) is locally given by 

Assume that the cost function is control-quadratic in the following way L = -(ul+ul). Prom ([SSD, §6 
we have 



iz^dt^ = S^ (Ci)i = l, (Ci)2 = 0, (C2)^ = 0, (C2)2 = l, 



dpf ^ {EAy = v' 



A' 



d,?^ — ^ (Gi)i + (G2)f = -A8(nicos6l + n2sin6'), 
dg2 — ^ (Gi)i + (G2)i = -A7(uicos6' + n2sin6'), 
dq' (Gi)H(G2)i = Ai + A7, 

dg^ ^ (Gi)i + (G2)l = A2 - As, 

d?' ^ (Gi)| + (G2)i = A3, 

d?' (Gi)i + (G2)i = A4, 

dt;^ ^ Ai + = 0, X2+pI = 0, A3 + = 0, A4 + pI = 0, A5 + pi = 0, 

-A5+pi = 0, Xq+pI = 0, -X^+pI = o, _A7+p2 = o, -A8+pi = 0, 
pi = 0, pI = o, 
dui — > Xyq"^ cos9 + Xgq^ cos9 — ui = 0, 
du2 — > Xjq^ sin9 + Xsq^ sm9 — U2 = 0- 
Hence all the controls and Lagrange multipliers are determined: 

Ai = -p\, X2 = -pI, A3 = -pf, A4 = -pi, A5 = -pI = pI, Ae = -pI = pI, A7 = pj, 
Xs = pI, ui = p^q"^ cos 6 + p^q^ cos 6 , U2 = p'^q'^ sin 6 + p^q^ sin 9. 
The cost function can be written as follows: 

l^=\{ul + ul) = {plq'+plq'f. 
Since P3 = and pi = 0, we have ZAipl) = Za{p\) = for ^ = 1,2. Then, {Gi)l = {G2)l = 
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(^"2)4 = = 0. Having this in mind, we have 





-(G2)? - 






-(G2)i - 


piipii^ +p^Q^ 


(G2)i = 


-Pi +Ph 




(G2)I = 


-vl-ph 




(Gi)^ = 


-{G2)l - 


pI, 


(Gi)^ = 


-{G2)l - 


pI 



(44) 



Note that the controls satisfy the following relationship ui sin 6 = U2 cos 9. If we impose the tangency 
condition, we have 

Zi{uism9 — U2Cos9) = (1^1)1 sin ^ — (1)1)2 cos ^ = 0, 

^2(^1 sin — ti2 cos 0) = ni cos + ti2 sin + (1)2)1 sin — (1)2)2 cos = 0. 

Thus, {Di)i = cos 9, (1)1)2 = sine and (1)2)1 sin 6* - (1)2)2 cos = -pW - pW ■ 
By conditions ([I6D, ([IT]) we have that {GA)f = for ^ / 5. 
After imposing the tangency conditions in (|4ip we obtain 

Za{^^) = {Fa)\ -4 = 0, Za{^^) = {Fa)1 -v\ = 0, 

Za{^') = {Fa)\ -v\ = ^, Za{^'') = {Fa)1 -4 = 0, 

Za{^'') = {Fa)1 - {Fa)\ = 0, Za(^6) = {Fa)\ - {Fa)\ = 0, 

Za{^'^) = -v\ - {Fa)2 + {F>A)iq'^ COS 9 - 5^uiq^ sin 9 + v\ui cos 9 + (DA)2g^ sin 9 + 5^U2q^ cos 9 

+v\u2 sin^ = 0, 

Za{^^) = v\- {Fa)2 + {DA)iq^cos9 - 6^uiq^ shi9 + v\uicos9 + {DA)2q^sm9 + 6^U2q^cos9 

+v\u2 s'm9 = 0, 
Thus, 



+ (^1 + (-Ca)i9^ cos 61 — SnUiq^ sinfl + w^ui 008 6* + (1)^)2'?^ sin 6* + (5^U2(7^ cos6' + u^U2 sin 6*) — — r- 
+ + (-Ca)i9^ COS0 — 82 "^iq^ sin 6* + d^ui cos6' + {Da)2Q^ sin^^ + (5^U2g^ cos6' + sin 



+ {FA)l^ + {FA)\^,+{GA)t^, 

ovl ov{ dpf 

where =t,t^ = 9, (Di)i = cos 6*, (Di)2 = sin0 and (1)2)1 sin 6* - (1)2)2 cos 6* = -Pqq'^ - pIq^ and 
also equations (|44p are satisfied. The optimal sections are integral sections of Z = (Zi, . . . , Z/^). 

6 Future work 

After this first geometric approach to optimal control problems governed by partial differential 
equations, it remains open to find the way to successfully extend any control system regardless of 
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the nature of the cost function. The main difficulty is to obtain a compatible system of partial 
differential equations after extending the original control system. 

In this paper we have not mentioned the different kind of extremals for optimal control problems. 
There exist the so-called abnormal extremals which are characterized at first without considering 
the cost function. As shown in [3], the constraint algorithm in the sense of Gotay-Nester-Hinds 
is useful to characterize the different kind of extremals in optimal control theory. Now, that the 
optimal control problems governed by partial differential equations have been understood in the k- 
symplectic framework, it seems that the application of the constraint algorithm for fc-presymplectic 
Hamiltonian systems [8] will characterize the extremals of those problems. 
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